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ABSTRACT 


The present day EHV/UHV transmission systems need a 
fast, reliable and efficient protective relaying system. 

In the past, electromechanical relays were being used. 
Due to their numerous disadvantages, they were replaced in 
the beginning by electronic relays and later by solid state 
relays. The advancements in IC technology have helped in the 
design of many compact relaying circuits. The advent of 
programmable equipment has revolutionarised the art of 
relaying . 

A relaying scheme for the protiection of a three-phase 
transmission line requires 6 measuring units to cater to all 
possible types of faults. 

In an attempt to reduce the number of measuring units 
polyphase distance relays were proposed. 

In the present thesis, a polyphase distance relaying 
scheme for the protection of transmission lines is proposed. 
The block diagram and the circuit configurations have been 
given. The proposed static polyphase distance relaying 
scheme has been fabricated and tested. 

Digital protection schemes for three phase and six phase 
transmission lines are also proposed. The proposed digital 
protection schemes have been tested on sample power system 
networks. 
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Chapter One 


INTRODUCTION 


1.1 MOTIVATION 

The art and science of protective relaying in power 
systems has assumed a greater importance all over the world 
with the tremendous expansion of interconnected power 
system networks. Thgire has been a demand for more reliable, 
efficient and fast acting systems of protection. This has 
led to the development of static relays and also digital 
protection schemes, which offer, to some extent a solution 
to the above mentioned problems. 

In the continued development of improved relaying 
schemes, the major targets have been increased speed and 
improved reliability. Of these, speed offers the greater 
economic contribution to the system. 

Studies have shown that a given system can transmit 
larger amount of power without loss of synchronism when 
short circuits are cleared in a .shorter time. Vi i th today’s 
static relays and fast acting circuit breakers, it is 
possible to remove a three phase faults on a transmission 
line in a maximum of 2 cycles and thus, maximise the systems 
load transfer capability. Conversely, the system may be 
designed more economically possibly with fewer lines, or 
switching stations when faster relays and circuit breakers 
are used. 
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Higher speed clearing also reduces burning of line 
conductors which may otherwise result in eventual conductor 
failure by subsequent corrossion with resulting loss of 
service. A second effect of pitted conductors is to 
increase corona and radio noise, which adversely 
affect the public communication systems. 

The large number of static relays in service since the 
year 1960 to the present have confirmed the high degree of 
reliability which can be achieved through conservative 
design of solid state circuits and complete protection 
against electrical surges induced from the neighbouring 
conductors. 

In the past O.C. relays were being used for the pro- 
tection of transmission lines. But owing to its inherent 
limitations such as, change of balance point with the type 
of fault, changes in generation, they were replaced by 
distance relays. Later on many high speed distance relaying 
schemes with different threshold characteristics were 
proposed. 

A distance relaying scheme for the protection of a 3-phase 
transmission lino requires the use of six measuring units. In 
an attempt to reduce the number of measuring units, polyphase 
distance relaying schemes were proposed. A polyphase distance 
relay is a single unit which is capable of protecting a 
transmission line against all the types of faults (whi*h are 
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ten in number for a 3-phase transmission line). These 
essentially being static devices, have low VA burdens. 

They have the additional feature of being inherently 
immune to power swings. 

Thus, the polyphase- distance relaying schemes, in 
addition to minimising the number of measuring units, also 
reduce the circuit complexity, since additional circuitry 
for inhibiting the relay operation during power swings 
is avoided. Now, with the advent of programmable equipment, 
many fast and efficient relaying algorithms have been 
developed. 

1.2 PROBLEM DEFINITION 

Development work on static relays started during the 
year 1945, initially utilizing thermionic tubes. With tho 
advent of transistor and other solid state components 
during the yeais 1950-60, the static relay development 
incorporated these devices. In the initial stages, 
static circuits were mere substitutes of their electro- 
mechanical counterparts. In the later stages, new 
principles have been evolved on hybrid and other unconven- 
tional devices and characteristics. 

Static circuits have been found to be well suited to 
the development of polyphase relays giving directional, 
distance and differential characteristics. 
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Later, phase and amplitude comparators gave place to 
instantaneous sampling comparators. Digital and analog 
computers entered tho field of system protection and 
revolutionarised tho concept of relaying. 

In the present thesis, a polyphase distance relay 
based on phase comparison principle is proposed. The 
relaying scheme has been designed and fabricated. 

Schemes for the protection of 3-phase and 6-phase 
transmission lines using a digital computer have also been 
proposed. The digital computer programmes have been 
tested for all possible fault conditions, on a sample 
system on the DEC- 1090 system at IIT Kanpur. 

1.3 HISTORICAL TREND 

The origin of polyphase distance relays dates back 
to the years 1954 to 1959, with the introduction of the 
induction cup type of relay for all phase to phase faults 
and then extended to tho static form for tho phase to 
ground faults. 

Sonemann et al proposed a polyphase distance relaying 
scheme for the protection of transmission lines [l], the 
main drawback of this relay was its inability *to provide 
protection against single line to ground faults.,. 
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As a complement to it, Rockfeller [2] developeda ground 
distance relay based on amplitude comparison of fault 
point sequence voltages. Another principle of polyphase 
distance relay by phase sequence detection of compensated 
phase voltages, shoved narrow operating zone on the complex 
impedance plane, specially during faults involving ground. 

A polyphase distance relaying scheme based on multi- 
input selective phase comparison was proposed by Choudhuri 
at al [4] in 1973. 

Although the present polyphase distance relaying 
schemes are successful in application, they suffer from a 
number of demerits such as, lack of flexibility, duplication 
of specification effort, lack of adaptability to system 
changing conditions and exclusion of self testing (in off 
line and on line) facilities due to complexity and cost. 

These disadvantages have resulted in a trend towards 
the use of programmable equipment in place of hard wired 
system. There has been enough work in this area and a 
number of fast and efficient algorithms have been developed. 

1.4 SALIENT FEATURES OF THE WORK 

The proposed static polyphase distance relaying scheme 
discussed in this thesis is capable of detecting any type 
of fault and producing a tripping signal in about half a 
cycle. It utilizes four 90° transient free phase comparators 
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based on block average phase comparison principle. 

Out of the four comparators, one is meant for catering 
to three phase faults, and the remaining three respond 
correctly for all single line to ground and all double 
phase faults with or without ground. The outputs from 
the four units are fed to a logic ’OR* gate which initiates 
tripping. 

It has been demonstrated in the above scheme that a 
ground distance relay utilising the shift between the 
phase to neutral voltages at the relay point and the same 
voltage properly compensated by line currents flowing 
through the replica impedances, can be made to respond 
correctly to both ground and phase faults. 

In the proposed digital relaying scheme, a digital 
computer program capable of detecting all possible types 
of faults and producing a tripping signal in less than 
a i- a cycle has been developed. A quadrilateral threshold 
characteristic which has been realised in the above 
scheme, has the following salient features. ' 

1. It is in operative during power swings. 

2. It is inherently directional. 

3. It is in operative at higher power factors. 
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1.5 CHAPTERWISE DESCRIPTION 

Now, a chapterwise description of the present work 
is given briefly. 

In Chapter 2, the general principle of polyphase 
distance relays is discussed* It is followed by the theory 
of the proposed relaying scheme and a description of the 
experimental circuit. The chapter concludes with the 
results and discussions. 

Chapter 3 deals with the general principles of digital 
protection. In addition to this, the merits of digital 
protection and different algorithms for real time calcula- 
tions of impedance are discussed. 

Chapter 4 deals with the theory and mathematical formu- 
lation including the algorithm for the digital simulation 
of the proposed three phase relaying scheme. 

A relaying scheme for the protection of a six-phase 
transmission line is also given in this chapter. 

Chapter 5 deals with the findings of the present work 
and scope for future exploration. 
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Chapter Two 

PROPOSED STATIC POLYPHASE DISTANCE 
RELAY 

2 . 1 INTRODUCTION 

A polyphase distance relay is preferred to single 
phase distance relays for the protection of transmission 
lines. The use of a polyphase distance relay reduces 
considerably the number of relay measuring units required. 
Further, the polyphase relays are not affected, in general, 
by power swings, momentary overloads and faults occuring 
behind the relay location. 

Switched relays were the earliest type of polyphase 
relays proposed. The distance measuring units in these 
relays are generally connected for phase faults (i.e. 

Delta voltages and Delta currents) and are switched 
to star connections when a single line to ground fault 
occurs. Such relaying scheme could not become popular due 
to several drawbacks such as complex switching circuits, 
over reach and under reach. 

Later on, polyphase distance relays using the poten- 
tials at the relay location and compensating for IZ drop 
in the line between the relay and the fault, were proposed. 
The nearest approach to a practical solution was a phase 
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phase comparator, compensated for line drops and zero 
sequence components to give the phase to neutral potentials 
at the fault which generate a zero sine product for a fault 
at the balance point. But this scheme had the disadvn 
adisadvantage of under reaching for a balanced 3-phase fault 
and thus required a second relay to clear such faults. 

Amplitude comperators seemed to offer a solution to the 
above problem. They could use either the phasor quantities 
in star or sequence quantities but neither of the two 
could give much economy over single phase relays and the 
use of sequence quantities introduced some loss in accuracy 
because of sequence filters. Also, the advantages of using 
three comparators instead of the six required for single 
phase relays is more or less cancelled by the cost introduced 
by the filters and complex circuitry. 

Later on, a method involving the determination of phase 
sequence of the compensated phase voltages was proposed. 

But it had the following disadvantages i 

a) the relay may not operate when there is a high arcing 
resistance near the balance point, 

b) the relay may maloperate due to DC offset present in the 
fault current, 

c) the operating time can be as high as 4 cycles for a fault 
near the balance point due to high DC offset current. 
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In the present thesis, a polyphase distance relay 
based on phase comparison principle of the compensated phase 
voltages is proposed. 

This relay does not have the short comings of the 
existing polyphase distance relays as mentioned above. 

2.2 THEORY OF THE PROPOSED POLYPHASE DISTANCE RELAYING SCHEME 

In the proposed scheme, a polyphase distance relay 
capable of detecting correctly all the single line to ground 
faults and all the double phase faults with or without ground, 
by the phase comparison of the line to neutral compensated 
fault point voltages has been realised. To cater for all the 
types of faults, four 90° transient free comparators based 
on block average phase comparison principle [4] has been 
used. Out of these four comparators, one is meant for 
catering to three phase faults and the remaining three 
respond for all single line to ground faults and double 
phase faults, with or without ground. 

The outputs from the four comparators are fed to a 
logic 'OR' gate which initiates tripping. 

Deri vation of the i nput quantiti es 

The line to neutral fault point voltage of any phase 'a* 
during any type of short circuit and assuming zero fault 
impedance is given by : 
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V fa V ra ” 

whore 

Vf a = compensated phase voltage of phase 'a 1 
V ra = volta! 3 e of P.T. secondary 
= line impedance, 

I a = line current in phase ’ a f * 

V fa = V ra -I al Z Ll _I a2 Z L2 _I ao Z Lo 
= V ra-< I a +KI ao> Z Ll 

where 

K = n-1, n = z Lo / 2^- and assuming = Z^ 2 for lines. 

Hence, the three compensated phase voltages at the relay 
location will be given by 


V x = V ra-<VKI ao )Z R - V xl+ V x2+ V x0 

V y = = a2v xi +aV xl +aV x2 +v 

V z = V rc -(X c+ KI a0 ) Zft . aV xl+ a 2 V x2+ V xo 


xo 


whore 

Z 


^ = Z r ,the replica impedance 

, 1 , V x2 and V xQ are the positive, negative and zero 


V 


sequence components of V . respectively. 

/C 

V xl - ha-hlV V x2 = V ra2 - Wr 
V xo = V rao - I ao 
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Fig. 1 gives the circuit diagram for obtaining V , 

A 

V y ,V z etc. 

i222£iE£2i3l_Circuit 

For the construction of polyphase distance relay, 

Q 

90 phase comparators using integrated circuits have been 
taken up. 

Relevant literatures [5] show that there are three basic 
methods, namely, 

a) Block instantaneous comparison 

b) Pulse comparison 

c) Block average comparison, a development of(a). 

In the block instantaneous comparison, the duration of 
polarity coincidence determines the output. The tripping 
criterion is that the duration of the first coincidence 
should exceed a specified time, usually one quarter of the 
power frequency period. 

In the pulse comparison method, the polarity of one 
signal is measured during a short interval in one cycle 
of the second signal usually, but not necessarily, at the 
latter's peak. 

In the block coverage comparison method, the duration 
of polarity coincidence is measured on both -the half cycles 
of the Input signals and the average value is determined in 
an integrating circuit, a trip signal being produced if a 
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specified average value is maintained for more than a 
prescribed duration. 

In the proposed scheme, block average comparison 
method is used. The two input quantities of the 
type jV x , V etc. are compared in a coincidence circuit, 
in which the duration of polarity coincidence is measured 
on both the half cycles of the input signals, in order to 
avoid transient over reach which might occur in (a) and 

(b) because of the assymetry in one of the signals* 

\ 

The standard dutput pulses are positive when jV ,V 

x y 

etc. are of the same polarity and are negative when they 
are of opposite polarity. The output pulses from the 
coincidence circuit are. fed to a linear integrator, whose 
output increases linearly when the pulse is positive 
and falls at the same rate when polarity reverse as shown 
in Fig. 4. 

When the integrator output exceeds some preset value, 
the level detector sends a tripping signal to the final 
tripping circuit, but the level detector resets when the 
integrator output falls below the preset value. 

The block average principle is preferred, because, 
although under steady state, there is no basic difference 
between the three, but under dynamic conditions, comparators 
in the categories (a) and (b) are inherently more responsive 
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to system transients and other spurious bysignals by 
virtue of their instantaneous operation, and their 
measuring accuracy can be maintained without sacrificing 
their speed. 

The rise and fall rates are at the designer's disposal, 
so that, the critical phase angle may be set to any 
desired value. Both the level detector set and reset 
levels are critical relation so the total excursion 
limits of integrator linearity and also to the slope of the 
output. Considding first the setting, it should atleast 
exceed a value which would be reached after one quarter 
of the system periodic time. If this were not so, the 
output would switch at twice the system frequency, even 
if the displacement between the input signals was greater 
than the critical value. The difference between set and 
reset levels should also exceed this same value in order 
that cyclic switching does not occur for marginal phase 
displacements when the net rate of change of integrator is 
very small. Finally, the upper limit of linearity should 
not be excessive, otherwise, the reset time will be poor. 

If all those factors are taken into account, together with 
the problem of designing a trigger circuit to operate to a 
specified level, it is found that the optimum level detector 
setting is 2/3 of the the integrator excursion limit and 
reset l/3 of the same limit. The block average comparison 



15 

scheme is inherently transient free, with a minimum 
operating time of 10 ms. 

Testing. 

Due to the nonavailability of a polyphase relay test 
bench, • the relay was tested as a phase comparator by 
giving two signals, representing the compensated phase 
voltages and it was found to give a tripping pulse over 
a band of + 90°. 
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Chapter Three 

DIGITAL PROTECTION OF TRANSMISSION LINES 
3.1 INTRODUCTION 

The digital protection of a transmission line has 
attracted much attention from protection engineers in 
recent years. 

Growing interest in the use of microprocessors for 
protection functions is evident from the. number of technical 
papers decreasing the study of and experimentation with 
such protection techniques. 

The advantages of protective relaying using digital 
computers are : 

a) Flexibility : 

A single general purpose hardware base, can be used to 
obtain a variety of threshold charactei istics and other 
functions with the change of the- software programming only. 

b) Adaptive capability s 

The processor can be programmed to automatically change 
its behaviour depending upon external circumstances which 
change with time. 
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c) Data interface access ; 

A general purpose digital computer systems can always 
be equipped with input/ouput ports through which data 
and control commands are exchanged. 

d) Self checking ability ; 

Unlike conventional relays, a digital processor is by 
nature a dynamic device and most hardware failures are 
flogged as soon as they occur by a processor stop. During 
normal operation, when there is no fault, specific 
programs can bo executed which also test the processing 
hardware. 

3.2 GENERAL PRINCIPLES OF DIGITAL PROTECTION 

Although a number of arrangements are possible for 
application of a digital computer to relaying, a simplified 
hardware configuration is shown in Fig . 5 as applied for 
the protection of a EHV transmission line. 

An analog input subsystem accepts 3-phase ac quantities 
from power system through conventional CT’s and PT's. All 
of these quantities are sampled simultaneously at uniform 
intervals, converted to digital form and transferred to the 
digital processor. The processor stores organises and 
makes decisions based on the values of the samples. Its 
prime purpose is to output a breaker-tripping command when 
a fault occurs on the protected line. Other secondary 
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controls and data output may also be provided. Stored 
in the processor memory, is an elaborate programme, 
which implements the line protection sensing and logic. 

At the core of the relaying program, is an algorithm 
which operates on raw incoming data samples, to provide 
a meaningful indicator of fault presen as well as its 
location. 

3.3 REAL TIME CALCULATION OF IMPEDANCE 

If the computer is fed with the sampled values of 
voltage and current signals, derived from the line, it 
will be possible to compute the peak current, peak voltage, 
phase angle, and impedance in each cycle. It is necessary, 
for the computer to store the sample values atleast over a 
period of 1 cycle of the system frequency. With these 
values, the digital computer computes all the relevant 
information. With the arrival of the new sample, the 
first sampled values in the previous ensemble is deleted 
and the same computation is again made over the present 
sot of samples. This process is continuously repeated by 
the on-line computer. 

The calculation of impedance based on the peak values 
and the phase angles based on the zero crossings will be 
erroneous because of the presence of transients in line 
voltage and current immediately after a fault. A few 



19 


methods for the calculation of impedance and phase angle 
are discussed in the following section. 

3.3.1 Imp edan ce calc ulation based o n peak valu es 

Among the sample values of voltage and current over 

a period of one cycle, the positive and negative peak 

amplitudes are found and their average magnitude is taken 

to represent V or I • 
max max 

V - i ^peak positive! + l^peak negative! 
max 2 

T - ^peak positive^ + f^peak neqative l 
X max ' ^ 

Z = V /I 

max' max 

Similarly, the time differences of V peak positive and 
Vale positive (I positive ) and that between V peak negatlve 

and Vak negative (T negatlve> are com P uted for each c Y cle - 
The average of these gives a measure of the phase angle 

<p ha , = [ T PQsitive + T negativej w 

where w is the system angular frequency in rad/sec. 

These computation are repeated for each cycle of the 
faulted waveforms. 
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3.3,2 Iffgodanc e calcula tion usi n g the integral s of vo l tag e 

If the voltage and current signals are assumed- to be 
sinusoidal 

v = V sin wt 
m 

i = I sin (wt - <p) 

then the integrals over a period of one cycle of rectified 
waveforms are given by 

V dt = 4V 
m 

i dt = 41 

m 

The effect of higher frequencies present in the 
transient waveforms will be almost nullified due to inte- 
gration. Thus the impedance measured in this way will be 
more accurate than the one measured from individual samples. 


2-n+a/w 

/ 

a/w 


( 2n+ ~) 


a/w 

for any a. 



21 


3.3.3 Transmission_line_gro tection_based_on system parameters 

This assumes representation of a line by a set of 
differential equation. The most common modes of a trans- 
mission line is described by the following differential 
equations 


v = R i + L § (3.8) 

The above representation of a transmission line recognizes 
the D.C. offset as a valid part of the solution and there- 
fore no special features need be implemented to suppress 
the D.C. osoffset. 

A number of algorithms have been suggested to solvc- 
equation (3.8) numerically. In 1971, Mcinnes et al. [14] 
proposed an algorithm for this purpose. It proposed 
integration over two successive time intervals, so that a 
sufficient number of equations are obtained to solve 'R' 
and ’L'. Integral equations are solved numerically using 
the trapezoidal rule. 


The final expressions for R and L take the form 


L 


R 


h 


= 1 


( V k _ 1 +V k ) ( ) “ ( V k _ 1 +V k __2 ) ( x k -i k _ 1 

( x k .jjt x k ) ( x k ^ x k - 2 )-(i k _ ^+x k _2 ^ ^ 1 i c “ x^ ^ 

[ V k-l +v k-2^ [ x k-l +i k^~^ V k-1 +V k^ 1 k-l +i k-2 
( i k _-j+x k ) ( x k _ ^*"X k ^ 2 ) — ( i k __^"ti k _ q ) ( i k — i k __.^ 


) 

) 

) 

) 


(3.9) 

(3.10) 
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whore ’v’ and ’ i* arc instantaneous values of voltage and 
current, 'K' is the instant and 'h 1 , the time interval. 

However, it should be noted that there are several 
problems, associated with the characteristics of actual 
transmission lines which are not accounted in eqn. (3.8). 

This equation assumes perfectly transposed lines and 
neither the shunt capacitance nor the series compensation 
is considered. 

The fault resistance and effect of power flow on the 
line at the moment of a fault is also not considered. The 
data window is one cycle. 

In 1975, Ranjbar et al. [8] developed another technique 
which relates appropriate integration intervals of equation 
(3.S ) to the particular harmonics, that are selected for 
removal. The sampling rate is related as a multiple of 
the order of the harmonic to bo removed, which makes the 
procedure quite restricted by the sampling rate selection 
and accuracy. 

3.3.4 Transmission line p rotec ti on bas ed on fundamental 
frequency sig na ls 

This relies on the theory of orthogonal transforms [15]. 
The most widely used in the fourier transform theory which 
utilizes the set of sine and cosine functions as a sum of 
combinations of the functions from the derived orthogonal 
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set. Basic properties of the fourier transform can be 
used to extract any particular frequency component from 
the incoming signal. In the case of the continuous 
functions, numerical approximation is done to obtain 
a digital solution. 


V 

I 



Fig. 5 Simplified Hardware Configuration 
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3.3.5 j;°U£^g2LJ!5glygj: | S method 

In this method, the ensemble of samples over a 
period of one cycle is assumed to repeat periodically and 
fourier analysis is performed on the ensemble of samples. 
The amplitude and phase angle of the fundamental component 
is obtained as follows : 


a l = 00 + 2^2 cos 2X + ... + f m cos mx) 

t>l = X/ 2 (2.f^ sin X + 2f^ sin2X -+ •... 2f m _^ sin(m-l)X) 


where f Q ,f^ ... f^ are sampled values of input signal over 
a period of one cycle and X is the sampling interval given 
by 2inc/m. The factors cox X, sin X, . .., cos mx, sin(m-l)X 
are constants and can be calculated and stored in the memory 
apriori as weighting functions on the sample values. The 
fundamental quantity is given by 


f x ( t) 




sin(wt+tan 


-1 



If this computation is made for both voltage and current, 
the impedance and the phase angle can be evaluated. 


3, 3 J6 PREDICTIVE CALCULATION OF PEAK VALUES 

In this method, the peak values of the voltage and 
current waveforms are calculated from each sample. 

!f v = V m sin wt 
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. then 

V m = f v 2 + ( V '/w ) 2 

and 

x m = f^+ TiVw)^ 

whore v and i are the voltage and current samples and v* 
and i’ are their time derivatives. These derivatives are 
computed numerically. 

The phase angle is given by 

? = [ten * 1 fi - tan ' 1 $£ ] 


Out of the above mentioned methods, the fourier 
analysis' method yields the best result, but it requires 
atleast a full cycle of samples of line voltages and 
currents immediately after the fault. 

In the proposed relaying scheme, the predictive 
calculation of peak values is used since it requires just 
3 samples to determine the peak values of voltage and 
currents if the central difference theorem of differentia- 
tion is used and only two samples if the backward 
difference theorem is used. 
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Chapter Four 

PROPOSED DIGITAL PROTECTION SCHEMES 

4.1 INTRODUCTION 

In the proposed relaying scheme a three phase relay 
with quadrilateral threshold characteristics has been 
simulated. A digital computer can be used to implement 
this feature on a practical system. The relay might 
prove its utility for on-line protection of present day 
EHV systems. 

4.2 THEORY 

The principle underlying this scheme is the detection 
of disturbances, i.e. faults, classification of the type’ of 
fault (i.o. phase or ground) and generation of ausuitable 
single phase set of relaying quantities such as voltage and 
current which are subsequently used for impedance calcula- 
tion. 

Three line to ground voltages and three line currents 
are scanned sequentially by subroutine 'SAMP'. Computation 
of zero sequence current and voltage is interleaved with 
the scanning and the subroutine 'SAMP' stores these eight 
quantities , each time it is called for. 

A counter, such as COUR for Red; COUY for yellow and 
COUB for Blue is assigned to each phase voltage. 
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The voltage samples in each of the phases are compared 
with the corresponding sample in the previous cycle by COMP. 

If the values differ in excess of a tolerance TOLV, the 
counter for that particular phase is incremented. If any 
of the counters exceeds a value CSV, the system enters into 
the fault detection mode. On the other hand, if the difference 
between two voltage samples found by the above comparison is, 
less than TOLV, the counter for the corresponding phase is 
decremented, if it is not already zero. 

When atleast one of the counters, COUR, COUY, COUB 
has a value equal to 4, the system enters into fault detection 
subroutines . 

4.3 CLASSIFICATION OF TYPE OF FAULT 

This is done by subroutines FTYPE and SELCT. The aim of 
these subroutines is to determine which two of the following 

1. Red phase, 

2. Yellow phase, 

3. Blue phase, 

4. Ground* 

are involved In the fault and to derive equivalent single 
phase relaying quantity. 

A counter for zero sequence voltage CZSG is used to 
determine if a fault detected involves ground. The last five 
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zero sequence voltage samples are compared and for each value 
exceeding a specified tolerance TOLZ, the counter CZSG is 
incremented. If the value of CZSG is atleast 4, it indicates 
that ground is involved in the fault. Subroutine SELECT 
determines whether any of the red, yellow or blue phases with 
or without ground is involved in the fault on the basis of the 
value of these counters. 

If a counter has a value 4, then the phase it represents 
is considered to be involved in the fault. If the counter 
has a value 2 or 3, the result is undetermined. The main 
aim of SELECT is to determine two of the counters for subse- 
quent use, i.e. if a RYG fault occurs, it would be sufficient 
to treat it as RY, RG or YG fault but in any case only two of 
the red, yellow or blue phase or ground need be considered 
as being involved in the fault. The red, yellow or blue phase 
counters of the same value as the ground counter are to be 
selected in preference to it (i.e. ground), since subsequent 
impedance calculations are more accurate, i.e., it is preferable 
to calculate a RYG fault as an RY fault rather than RG or YG 
fault. In case of a single line to ground, fault, zero sequence 
compensated current is used for subsequent impedance calcula- 
tion. 
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4*4 MATHEMATICAL FORMULATION FOR IMPEDANCE CALCULATION 

In the present scheme, it is proposed tow rwork with 
delta quantities for phase faults and phase voltages and zero 
sequence compensated current for ground faults because point 
in this case does not shift due to different type of faults. 

The method of calculation of line impedance [5] involves 
the predictive calculation of peak current and peak voltage, 
the impedance being determined by division of peak voltage by 
peak current. A digital computer sampling a sinusoidal wave- 
form, can determine the peak values as they occur. 

However, it is necessary in the interest of time to 
determine the peak values before their occurrence, i.e. to 
predict the peak value of the waveform from any given sample. 


Let v be the voltage sample at any instant, and V ^ 
the peak voltage. 

Since it is not proposed to synchronise the sampling to 
the sinusoidal, the value of sin wt is also unknown. 


Lot v = V 


sin wt 


Differentiating (4.1) we get. 


v' 



cos wt 


(4.1) 


(4.2) 


The various numerical methods for differentiation are given 
in Appendix III . Central difference method, being the best 
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from the point of view of accuracy, is being used. 

From eqns. (4.2) we get, 

w“ = V P k cps wt (4.3) 

Squaring (4.1), (4.2) and adding, we ge t , 

(v’/w) 4 . (v ) 2 = v 2 k [cos 2 wt + sin 2 wt] 

V pk = (v'/w ) 2 + (v ) 2 ( 4 . 4 ) 

Similar ly peak value of sinusoidal current is, 

I pk = (i'/w ) 2 + i 2 ( 4 . 5 ) 

Dividing, equation (4.1) by (4.3), we get 
vw _ v pk sin wt 

V ' “ _ , 

V cos wt 

vw 

tan wt = 

v> 

wt = tan ” 1 [~] 

Thus point on cycle of voltage sample, 

V @ = arc tan(~) (4.6) 

Similarly, the point on cycle of a current sample, 

Iq = arc tan(jy) 


(4.7) 
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Impedance modulus by use of eqn. (4.4) and (4.5) is, 

,2 




Z 2 

& 


(4.8) 


Impedance phase angle is 


z e ~ I e 


'e 


= arc tan(ji) -arc tan(~f) 


wv • 
v 1 


(4.9) 


Resistance ’R f and reactance 'X f upto the fault point is 
calculated from Z and Zq. The flow chart for the proposed 
relaying scheme Is given in Figv 6. 

4.5 RESULTS 

The proposed three phase relaying scheme has been tested 
on a sample power system network [8], the details of which 
are given in Appendix II for different types of faults* 

In the present scheme, a quadrilateral characteristic 
has been realised. This scheme possesses the following 

features : 

a) It is immune to power swing 

b) It has inherent directional feature 

c) The characteristic can be varied by changing the slopes 
at the instant relay starts 

d) It possesses a very high degree of sensitivity, reliabi- 
lity and is fast in operation. 
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The result of the test is entered in the Table 1. 

It has been found that the predicted value of the peak 
voltage and current for all types of fufaults, is accurate 
to within -0.01818 to 0.3083 percent and 0.4386 percent 
respectively. The Z modulus is accurate within 
-0.45 to 0.256 percent and the argument of Z to within 
0.328 percent. 

The time required for relay operation is determined by 
four considerations as shown. 

Firstly, four samples are required for detection of 
occurrence of a disturbance. Secondly, seven samples are 
required for impedance calculations, these may include 
those required for detection of fault. Thirdly, the 
impedance calculation itself will take an additional time, 
and fourthly, the calculation of R and X after calculation of 
impedance, requires a finite time. 

The proposed relaying scheme which is a compact 
unit for the protection for a 3-phase transmission line 
is also capable of detecting the location of any type of 

fault on the line. 

The operating time of the relay was found to be 6 ms 
which is less than £ a cycle. 
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4.6 PROTECTION OF SIX PHASE TRANSMISSION LINES 
4.6.1 Introduction 

One of the most important items in applying the six- 
phase alternative in transmission planning is the design 
of an adequate protective scheme. 

In a six-phase system, the number of fault types is 
much greater than that of a 3-phase system-. The fact that 
there are six phases, each subjected to a different voltage, 
and a neutral in a six-gr&ase system, makes the number of 
fault types much greater than for the three-phase system. 
Out of the 120 possible combinations in a six-phase system, 
there are 23 combinations with distict fault levels and 
phase interconnections. 

Consider line to line faults that do not involve ground. 
The current magnitude and circuit assymmetry for a fault on 
phases a and b is the same as that on b and c and so on. 
Therefore, the 15 two phase fault combinations reduce to 
3 significant combinations. These are \ 

1) faults between phases 60° apart, e.g. a to b 

2) faults between phases 120° apart, e.g. a to c, 

3) faults between phases 180° apart, e.g. a to d. 

Similarly, the analysis for all faults can be confined to 
the 23 significant combinations. In the case of a 3-phase 
system, there are only 5 significant fault types. 
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4.6.2 Pr oposed prote ctio n s c heme for a 6 -phase trans m issi on 
line 

The proposed 3-phase digital protection scheme, in the 
previous section is extended for a 6-phase system* the only 
difference being in the number of counters, (seven counters 
C0U1, CCU2, C0U3, COU', C0U5, C0U6,for the six phases and 
COUG for the zero sequence voltage) and the type of 
characteristic realized (plane impedance characteristic in 
the present case). Explanations given in the previous 
sections hold good for the present scheme also. 

4.6.3 Classifica tion of t he t yp e of fault 

This is done by subroutines FTYPE and SELCT, The aim 
of these subroutines is to determine which two of the phases 
or ground are involved in a fault and to derive equivalent 
single phase relaying quantity. 

If a counter has a value 4, then the phase it 
represents is considered to be involved in the fault. If the 
counter has a value 2 or 3, the result is undetermined. The 
main aim of SELCT, is to determine two of the counters for 
subsequent use, i.e., if a 1-2-3-G fault occurs, it would be 
sufficient to treat it as 1-2, 2-3, 3-1, 1-G, 2-G or 3-G 

but in any case, only two of the six phases or ground need 
be considered as being involved in the fault. The six- 
phase counters are to be selected in preference to the 
zero sequence counter, since subsequent impedance calculations 
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are more accurate, i.e. it is preferable to calculate 
a 1-2-3— 4-G fault as a 1-2, 2-3, 3-4, 4-1, etc. fault 
rather than 1-G, 2-G, fault. In case of a single line to 
ground fault, tzero sequence compensated current is used 
for subsequent impedance calculation. 

4.6.4 Resul ts 

The proposed six phase relaying scheme has been tested 
on a sample system network, the details of which are given 
in Appendix II for different types of faults. The, results 
of the test are entered in Table 2. During the test, it, was 
observed that all the faults were detected satisfactorily 
i.e. in no case did the programme falsely indicate the 
involvement of the six phases and ground in a fault. 

It has been found from the test that for all the types 
of faults that the magnitudes of peak voltage and current are 
accurate to within 0.19506 percent and 0.17832 percent 
respectively. It has also been found that the impedance 
magnitude is accurate to within + 0.087 percent and the 
argument of the impedance to within -2.41 percent to 
0.3695 percent. 

The proposed relaying scheme, which has been tested 
on the DEC- 1090 computer system at IIT/Kanpur is a compact 
unit for the protection of a 6-phase transmission line. 

Tho operating time of the relay was found to be 6 ms 

which is less than % a cycle. 
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CHAPTER 5 

CONCLUSIONS 

The present work proposes protection schemes for 
both three phase and six phase transmission lines. 

The proposed polyphase distance relay offers 
protection against all types of faults likely to occur on 
a three phase transmission line. 

By the use of transient free 90° phase comparator 
based on block average phase comparison scheme and with 
simple circuits using ICs, a polyphase distance relay 
has been developed. The proposed static relay has been 
designed, fabricated and tested. 

The relay has the following advantages : 

1. It is not effected by arcing resistance 

2. Since, in the 90° phase comparators, the duration 

of polarity coincidence is measured on both the half 
cycles of input signals, the over reach in the relay 
is minimised, which may otherwise occur due to d.c. 
offset in the current waveform. 

3. The circuit is simple and devoid of any complex 
switching elements. 

The operating time of the relay has been found to be 
about a half a cycle. 
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The proposed polyphase digital protection scheme, 
in addition to responding correctly for all the types of 
faults, also possesses the most desired threshold chara- 
cteristics namely the quadrilateral characteristics. It is 
a flexible scheme which upon detection of a fault selects the 
phases involved and chooses the appropriate set of relaying 
currents and voltages for subsequent impedance calculations. 

The digital computer program computes the impedance to 
a reasonable accuracy, this is in spite of considering only 
the first term in the expression for numerical differentia- 
tion. Even better accuracies can be realized by considering 
a few more terms in the above expression. 

The proposed digital computer program has been tested 
in the DEC 1090 system, which showed a minimum operating 
time of about 6 ms which is less than \ a cycle for all 
types of faults. 

The six-phase protective relaying scheme developed 
is capable of detecting all possible types of faults on a 
six-phase transmission line. Even in the present scheme 
the impedance computation was within reasonable accuracy. 

The proposed scheme demonstrates a method for real time 
calculation of impedance for a six-phase transmission line. 
Upon computation of impedance, any complicated threshold 
characteristics can be easily realised as in the three-phase 


scheme 
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The operating time of the relay was also found to be 
less than y a cycle. 

Scop e for further work in_this_area 

The digital computer based relaying programme for the 
protection of three-phase generators and bus bars can be 
developed. 

A microprocessor based scheme for the protection of 
three-phase to six-phase transformers and buses can be 

developed. 

The computer application for on-line protection may 
lead to an overall computer control of power system networks 
in which a central computer would perform all the functions 
such as monitoring, control and protection. 



Table 1 (Three phase system) 
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Table 2 (Six phase system) 
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appendix I 


Transmission line data [8] which is used fop tosting ^ 
proposed protection scheme (3-phase). 

Base = 230 KV, 1000 MVA 
Length of line = 300 miles. 


Line parameters : 

Positive sequence reactance = j0.123xl 0 - 2 pu/mile 
Positive sequence resistance= 0.223XKT 3 pu/mile 
Zero sequence resistance = 1.226X10' 3 pu/mile 
Zero sequence reactance = J0.32xl0- 3 pu/mile 
Positive sequence capacitive reactance = tJ4 20 pu/mile 
Zero sequence capacitive reactance . -J714 pu/mile. 
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APPENDIX#' 


Transmission line data which is used for testing 
the proposed six-phase protection scheme* 

Base 100 MVA/phase 

Transmission line data 


Positive sequence impedance 
Zero sequence impedance 


= 0,0024 + jO.0116 p.u. 
= 0.0141 + j0.090l p.u. 



APPENDIX Ilf 


DIFFERENTIATION FORMULAE 

The basic central difference expression for the 
derivative [2] § 

hy£ = (p6- 6 3 + ^ 6 5 )y k (AII.l) 

where ^ jj, and 6 are the standard notations for the 
operations of backward differing, central differencing 
and averaging respectively and h is step size, and y stands 
for v or i. 

Using the first term only 

h y k = ^ ( y k+l “ y k-l^ (All. 2) 

and with a second-term 

hy k = - h (y k+2 + I y k+l * § Y k-1 + h y k-2 } (A 11 - 3 ) 

For backward differences 

hy£, = ($? - 2 - . . .)y k (All. .4) 

Using the first term only, 

hy k = y k- y k-i (AII - 5) 

and with a second term 

hy k = ^( y k- y k-2 } (AH. 6) 

Note that (All. 2) and (All. 3) for y k involve sample values 
at times later than t k , whereas (All. 5) and (All. 6) do not. 
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r 1 PHASE RELAY WITH QUA DR 1631 MSltl L CffAM2Cf»itTlT£&f 

{-***** + **# ******* **************************■$***********%******* 
r ************************************************************************* 

C wore 1 )=*COUNT£R FOR PHASE 1 

C NnF(2)=CnuNTER FOR PHASE 2 

r NtlF ( 3 ) ^COUNTER FDR PHASE 3 

C NOF(A)aCOUMTER FDR "©'SEQUENCE VOLTAGE 

C VIJZaPHASE VOLTAGE OR DELTA VOLTAGE DEPENDING UP ON THE TYPE OF 

c fault wether lg or any phase fault respectively, 

c TUZlsPHASE CURRENT OR DELTA CURRENT DEPENDING on the TYPE of faut.t. 

P$**!|^***4s*#$***f}fc#**##**$$;#******#;|t****#*********I*ff****jjt****#l|!*#*¥#i|S>|S**atE 

REAL K,IA,ITHETA,IUZ,VPEAK,lPEAK,l M P,IUZt ,TPEAK1,IT1 
INTEGER Z,W,'’,p,p,s, THETA, P i 

COMMO I V f 6 , 4 C J , AV C.6 , 40 ) , DV ( 40 ) , DI(40 ) , NOF (tO), 011(40) 

COMMON K(6»40),IAC6,40),VUZ(40),IUZC40),ZVC40) ,ZT(40),TUZ1C40) 
par-, f l«^■ ! I ,,, *?v,'D ^; 'VTCE^• # DSK # ,FILE• # LG,CDR , ■} 

; . ' ;v r ;;DO 222 J»l,36 

READ (20,223), C V a , J) , AV U , J) , I A (I , J) , Ial , 3 ) 

FORMAT (//9CF10.7,3X)) 

"Mi' CONTINUE 

ARA*,0S36lARB*,0fi6O?ARC».07?Xs,3fi9 

TY©E 6S 

65 FDP;.?AT t//3 oX, # 3 PHASE RELAY m/TTH DU AD CHAR TEST 

65 1 FOR A LLO FAULT"//) 

TY D t. 66 

66 FORMAT (30Y,3«JC '-")//) 

Wall 4.0 

406 D»J 1 1 = 1,3 

c Call ttmf(ni) 

1 N 0 F ( T ) a 0 

nqf( • 1 ) =o 
no 7 J = 1 ,36 

v~ i 

CALL S A M P ( P ) 

CALL COMP(P,TAF) 

TFCIAF.GE.1) GO ■ TO 3 

7 CONTINUE 

7 o-p+1 

?2 call sampc 0 ) jcall gompcp,iaf)icall ftypecp) ;tc«o 

Du 47 lal # 3;IFCN0Fm ,FQ,3‘) IC*IC+1 
TF(IC,LT,2,ANU.IC,NE,0) GO TO 3 

47 CONTINUE 

C TYPE 55 



*5 5 


FORMAT C///20X,'ALI VAT UEs l' ! PER tvtjv/) 

PsP+ 1 

FALL SAMP(F) 

T.i, FTYPEt P) 

CALL SFL^TCP) 
p = p + i 

CALL SAMP(P) 

CALL FTYPE(P) 

CALL SELCT(P) 

P = P+ 1 * CALL SAMP(P) JCALL FTYPE(P) |CALL SELCT(P) 

Hal ,/l80O. 

CALL TJMECN1) 

r>vcp-i)»cvuzcP)-vuzcp-?})/(2,*H) 

r>ICP*l)»(I!JZfP>~IUZ£P-2))/(2, *H) 
DIl(P-.U*£TUZl{P)-IUZl(P-2))/(2.*H) 
VPEAK»S0PT((VUZCP-1)**2)+DV(P-1 )**?/W**2) 

IPEAK»SQRTC£IUZC p -l ) **? ) +DI £P«1 ) **2/W**2) 
IPMiUi*«QRf£aOZt£PM)**2)+DlifP-n**2/W**2) 

T.F £CARS(VPEAK)).GE.1 ,414) GO TO 406 
IPC £ ABS £ IPEAK) ) , LE, 1 ) GO TO 406 

wv*w*vozcp-n;wiaw*iuz(p-n 

V THETA«ATAN20*V,nV£P-n ) 

TTHETA»ATAW2£WI,DICP-1)) 

TT1=AT»N?( 'T*TU7.1CP-n,nil (P-l )) 

VTaV’»’HPTA*t80,/3 f 1416 

n' r = TT'-»ETA*L«0/3. Hlfi 

7 T !f fi T A s V T HE T A «• T T 1 
Ti-|P = VWFA,F/TPE:AK.1 

PHT = 7T«KTA*l«-y,/3,t4l 6 
TFfPMl.LT.0,0) GO TO 406 
Nfl =T ( !D*ro, >Ht) ;Yt*IMp*STND(PHT) 

AKlsX/ARA; AKOaT/UROARB) 

T E C AR 1 , (jT , 0 , 0 ) GH TO 406 
T F ( X 1 , OT , , 6 ) GO TO 406 

TFCAPl.LF, AKR.A'in.Xl.LE, CAKl*ARl)) TYPE 100 
TFCAP1 ,GT. ARR.AnD.X 1.GE. (AK2*CARl-ARBn ) TYPE 100 

CAT.L TTME(M2) 

TYPE 45 

45 FORMAT C///20X,'PEAK VHlTAGE # , 5X , 'PEAK CTJRRENT', 

1 5X,'TMPE0ENCE',5X, 'PHASE A«GIjE C DEGREES )' , 

i 5 X » ' vqTiT amgt.e',5x, 'Current angle'//) 

TYPE 46, V PEAK , 1PE , I HP , PHI , VT , A XT 

46 FORMAT (?0 Y , r 7.5*10X,F7„3,10X,F7.5,inx,F7.1 , 1 6X , F? , 3 , JOX , FR . 3/ / ) 


TfP IT l *s 


■*’5 C2-:;y, # C^un:Pt? P*«l v ,«x , 'COUNTER PK2',5Y,'CniP'TFR p«3 # , 

I 5X,'0 SEO COUNTER'//) 

'”Y ?, E 36, f('MO p (T6) ,16*1,4) ) 

16 FORMAT C?OX,3(T3,l5X),lXrI3//) 

TY°E 17, P 

17 FORMAT (//4qX,'N0 OF SAMPLES REAP *',2X,T3//) 

IF ( »HI „GE„ 1*0,0, AND, PHI, LT, 360.) 00 in 406 
TFCIMp.GT.l ) GO TO 406 

TYPE 100 

100 F ORMA T (/////, 40 X*'TR IP CKT BREAKER') 

TFfWOFf 4).GE,4) f Y pe 44 

44 F0RMATC//40X, 'GROUND FAULT DETFCTEO'//) 

go to tot 

GO TP 101 

101 NR*N2»M1 f TYPE 133, NR 

131 FORM AT C//20X, 'OPERATING TIME OF THF RELAY s' , 1 X , T3 , 1 X , 'n>S '// ) 

STOP 

ENf) it- • , / /. >.■ • •* 

SUBROUTINE SAMP(P) 

THIS ROUTINE CALCULATES THE 'O'SEQUENCE VOLTAGES AND CURRENTS, 
C************************************************************************* 

PEAL K,IA,TTMETA,IUZ,VPEAK,IPEAK,IMP,IPEAK1 
INTEGER 7,,W,U,P,R,S,THETA,P1 

f o M m n i ; V(6,40),AV(6,40),PV(40),DIC40) ,NOF£10) ,0J1 (40) 

Ft jo,. KC6,4P) , !A(6,4o) ,VUZ(40) ,TUZ(40),ZVC40),ZI(40),IUZ1C40) 

? :>j ( P l s ( c v ( 1 , P ) + M ( 2 , P ) + V ( 1 , p ) ) ) / ( 3 , } 

71 nn = f (T a( 1 ,P) + TA(2,P)+lA(3,P) ) )/C3. ) 

RE T, IRN 

(7 D 

f * ************************************************************************ 

SUBROUTINE FTYRE(P) 

C TUTS ROUTINE SETS THE 'O'SEQUENCE COUNTER AND ALSO PREPARES ALL 

r the counters for routine 'sei-CT' 

r****** ***************** ************************************************** 

PEAL K,U,TTHETA,IUZ,VPEAK,IPEAK,IMP,IUZ1 ,IPEAK1 

INTEGER 7.,w,u,p,r,s,theta,pi 

COMMON V(6, 10) , AV(6,40) , DV ( 40 ) , Dl C 40) , NOE( 1 0) , 01 l ( 40) 

cnMfiOi-1 K(6,40),IA(6,4»),VUZ(40),IUZC40),ZV(40),ZI(40),TUZ1{40) 

integer t 

no 2 ? t»i , i 

TEf N n FCT) » G E • 4 ) CO TO 73 
K(T,P)«0,0 

GiJ TO 7 2 



r >2 '"’out tune 

P5asP»5 

PQ 21 Rxt>5,P 

TFf ARSCZV(R)-2V(R-1) ),LX, 0.02) GO tq 21 
A! U ^ ( 4 } a W 0 F C 4 ) + 1 

9 1 CONTINUE 

RETURN. 

Erin • 1 

C.*******'*****ni^*****t-M*-***'*t*-¥**.***-4*^* -Hi* ******** t* ************ *********** 

SUBROUTINE SRliCTCP) 

r SELECTS W-HICH Ttff) Op THE 3 PHASES oh GROUND &RF INVOLVED IN A FAULT 

C ANO TO DERIVE EQUIVALENT SINGLE PHASE RELAYING QUANTITY 

rt *************** ****: m******¥*************************************** 

REAL K, JA,ITHETA,inz,VPEAK,lPEAK,IMP f inzi ,TPEAK1 
COMPLEX Zl,Z2,ZO,PKl 

7/‘-,n,P r R,8,THPTA f Pl 

r D T '-‘0.' VC6,*0),AV(6,4u),DVC40),DTC4o),NOF(lO) r nil (40) 

COMMON KC6,40) ,IA(S,d(n ,VUZ(40>,TU7C40'> ,ZVC40) ,ZT(40> ,TU74C40) 

DIMENSION B(6 f 40),AC6,d0) 

Zl*C,0669 # .369);ZOa(.367S,.96) 

pKi*czo-zn/zi 

PK»CARS(PKt ) 

7,b0 ' ■■;■' ■■■■'' : . 

nn 2R Usi , 3 

°7 T F f K f !j , p 1 . N E . 0 . 0 ) GO TO 98 
GO TO 95 
78 7=9*1 

AC%P)=KCU,P1 

Pt7,o)=Lncn,P} 

95 ooWTTNUE , 

TFfZ.GE.7) GO T 0 100 

IFfZ.E0,1 .AHD.NOF(4),GT.O) GO TQ 60 
GO p | 

60 VU7( a )=Af 1,P) 

TU7.(P)»Bf l,P) 

To7 t CpIsPCI ,P) + ( D KIZT(P}) 

GO TO 61 

too VUZIB)*AM,P)-M-2,P) 

T |7,(P)aBCi ,P1-R(7,P) 

TU7.1 CPlaTUZ(P) 

61 RETURN 

Eon 



subroutine eohP(p,TAF) 

C noRS THE CYCLE by ryr,iw, COMPARISON OP VOLTAGE SAMPLER 

^*+**********♦***♦*****************♦#♦**************♦***#************♦#**4 
PEAL K,IA, TTHETA,IUZ f V»EAK,IPEAK,IMP f IUZ1 # IPEAK1,IT1 
t *TEGER 7,W,!J,P,R, s, THETA, PI 

COMMON VC6,40) ,AV(6,40),DVC40),DIC4Q) ,N0F(10),DI1 (40) 

COMMON K(6,40) r IA(6,40),VU7(40) f IU7,C40),ZVC40),ZT(40),TU7if40) 

IAF = 0 

o<J 20 Mat ,3 

IFC (ABS( AVCM,P)-Vf¥,P)),LT,0„05)) GO TO 10 
Ni.)F (, M ) *nOF ( M ) 1 1 ? GO TO ?0 

19 TF(wnF(M),NE.O) NQF ( M ) ayOf' ( M 5 "1 fclV" 

7 0 CONTINUE 

Du 2 1*1 f 3 

IF (N0P(T).GE,4) IAF*IAF+1 

2 CONTINUE 

RETURN LEND 



